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Abstract 

Aedes aegypti is the major vector of yellow and dengue fevers. After 10 generations of adult selection, an A. aegypti strain 
(SP) developed 1650-fold resistance to permethrin, which is one of the most widely used pyrethroid insecticides for 
mosquito control. SP larvae also developed 8790-fold resistance following selection of the adults. Prior to the selections, the 
frequencies of VI 01 6G and F1534C mutations in domains II and III, respectively, of voltage-sensitive sodium channel (Vssc, 
the target site of pyrethroid insecticide) were 0.44 and 0.56, respectively. In contrast, only G1016 alleles were present after 
two permethrin selections, indicating that G1016 can more contribute to the insensitivity of Vssc than C1534. In vivo 
metabolism studies showed that the SP strain excreted permethrin metabolites more rapidly than a susceptible SMK strain. 
Pretreatment with piperonyl butoxide caused strong inhibition of excretion of permethrin metabolites, suggesting that 
cytochrome P450 monooxygenases (P450s) play an important role in resistance development. In vitro metabolism studies 
also indicated an association of P450s with resistance. Microarray analysis showed that multiple P450 genes were over 
expressed during the larval and adult stages in the SP strain. Following quantitative real time PCR, we focused on two P450 
isoforms, CYP9M6 and CYP6BB2. Transcription levels of these P450s were well correlated with the rate of permethrin 
excretion and they were certainly capable of detoxifying permethrin to 4'-HO-permethrin. Over expression of CYP9M6 was 
partially due to gene amplification. There was no significant difference in the rate of permethrin reduction from cuticle 
between SP and SMK strains. 



Citation: Kasai S, Komagata O, Itokawa K, Shono T, Ng LC, et al. (2014) Mechanisms of Pyrethroid Resistance in the Dengue Mosquito Vector, Aedes aegypti: 
Target Site Insensitivity, Penetration, and Metabolism. PLoS Negl Trap Dis 8(6): e2948. doi:10.1371/journal.pntd.0002948 

Editor: Pattamaporn Kittayapong, Mahidol University, Thailand 

Received November 25, 2013; Accepted May 3, 2014; Published June 19, 2014 

Copyright: © 2014 Kasai et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This work was supported by Grants-in-Aid for Young Scientists (B) (23590501) and Grants-in-Aid for Scientific Research (C) (25460525), and was 
supported partially by grants for Research on Emerging and Re-emerging Infectious Diseases from the Ministry of Health, Labor and Welfare, Japan (H21-shinkou- 
ippan-005 and H24-shinkou-ippan-007). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the 
manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 
* E-mail: kasacin@nih.go.jp 



Introduction 

The yellow fever mosquito Aedes aegypti inhabits tropical and 
subtropical regions worldwide and is the major vector of dengue 
fever (DF) and yellow fever. DF is a rapidly growing health issue, 
with the average annual number of cases being approximately 1 00 
million and a 30-fold increase in the past 50 years [1,2]. The 
disease is endemic in at least 1 1 2 countries, especially in south and 
Southeast Asia, and an estimated 2.5 billion people are currently 
living in risk areas [1]. At present, DF causes more illness and 
death than any other arbovirus disease in humans [3-5]. 

Successful population control of vector insects is the key to 
prevent transmission and epidemics of infectious diseases. This 
strategy relies heavily on insecticides, particularly pyrethroid, a 
popular class of insecticides with high and rapid toxic activity 
toward insects and low toxicity to mammals [6]. Pyrethroids are 
used to control and/or prevent adult mosquitoes by ultra-low 
volume sprays, thermal fogging, pyrethroid-impregnated nets, etc. 
However, many dengue endemic areas are now facing the 
problem of pyrethroid resistance due to frequent and intensive 
use of these chemicals [7] . Resistance of A. aegypti to pyrethroids 



has been reported from various countries [8] . Understanding the 
level and mechanism of resistance to insecticides is essential for 
developing appropriate vector control measures. 

A. aegypti is present in most residential areas of Singapore. 
Despite a well-established national vector control program that 
includes community engagement, law enforcement, and inter- 
sectional coordination, Singapore continues to face the risk of DF 
resurgence [9-1 1]. During an outbreak in 2005, more than 14,000 
DF cases were reported [12]. Since 2005, surveillance for dengue 
control has been based on four pillars: (1) case surveillance through 
mandatory notification of dengue cases to the Ministry of Health, 
by all medical practitioners; (2) vector surveillance through 
premises checks by vector control officers from National Environ- 
ment Agency (NEA); (3) virus genotype surveillance at Environ- 
mental Health Institute of NEA; and (4) monitoring of other 
environmental parameters such as weather factors and population 
density. Clustering of cases and development of risk maps using 
the surveillance data allows prioritization of vector control 
operations. Though the program is largely based on source 
reduction and large scale fogging has not been conducted by the 
authorities since 2005, insecticides continue to be used by the 
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Author Summary 

Aedes aegypti inhabits tropical and subtropical regions 
worldwide and is the major vector of dengue and yellow 
fevers, and a secondary vector of chikungunya fever. 
Dengue fever is epidemic in more than 110 countries and 
causes up to 100 million infections annually. As there is no 
efficient vaccine or medicine currently available, vector 
control remains the primary solution for reducing the 
number of cases of this disease and relies heavily on the 
use of insecticides. Intensive and long-term use of 
insecticides has resulted in the worldwide emergence of 
mosquitoes with resistance to these chemicals. Here we 
newly identified two P450s, which have the ability to 
detoxify pyrethroid insecticide and are over produced in 
the resistant A. aegypti strain. Our study showed there 
were at least four P450 isozymes associated with resistance 
and target site insensitivity. These findings may lead to the 
development of more accurate monitoring systems in the 
field and also assist to identify new target sites for 
insecticides that are effective against resistant insects. 



authority for indoor misting in areas with dengue transmission, by 
private arrangement and household use of aerosol cans. Regular 
monitoring for insecticide has revealed that Singapore's A. aegypti 
larvae have developed resistance to synthetic pyrethroids, with 
resistance ratio to permethrin in the range of 29-47 times [13]. 

Generally, insecticide resistance in insects is caused by three 
major mechanisms: (1) reduced sensitivity of the target site, (2) 
reduced penetration of the insecticide due to altered cuticles, and 
(3) increased activity or level of detoxification enzyme(s). Insects 
develop resistance to insecticides by obtaining one or more of these 
mechanisms. 

Amino acid substitution in the voltage-sensitive sodium channel 
(Vssc), the target site of DDT and pyrethroids, is the most well- 
understood mechanism conferring resistance against pyrethroid 
insecticides [14,15]. This mechanism is called "knockdown 
resistance," with its inherited trait, "Mr" being first identified in 
the housefly as a leucine-phenylalanine substitution at position 
1014 (L1014F) of Vssc. Recently, many other mutations that cause 
reduced sensitivity to Vssc have been identified from many 
arthropod species and are called either "Mr-like factor" or 
"knockdown resistance." Knockdown resistance is unaffected by 
synergists that inhibit the activity of metabolic enzymes such as 
carboxyl esterases or cytochrome P450 monooxygenases (P450s). 
Several amino acid substitutions have been detected in A. aegypti 
Vssc: S989P, 1101 1M, 1101 IV, V1016G, V1016I, F1534C, and 
D1763Y [16-22]. Of these, only substitutions V1016G, V1016I, 
and F1534C have been shown to strongly correlate with 
pyrethroid resistance [17,23-27]. The dominance of these three 
mutations on pyrethroid sensitivity, however, is not well under- 
stood. 

Metabolism-mediated insecticide resistance is now considered a 
key mechanism in insects [7,28,29]. Three families of metabolic 
enzymes have been implicated in the metabolism of insecticides: 
esterases, glutathione transferases (GSTs), and P450s. The genome 
project of A. aegypti identified 26 esterases, 49 GSTs, and 160 
P450s [30]. Identification of the specific enzymes involved in 
insecticide resistance, however, has proven challenging. P450s has 
been shown to be the metabolic enzyme most strongly linked to 
the development of pyrethroid resistance [31-33]. However, due 
to the large number of P450 genes and the structural similarity 
among different isoforms, identification of the isoform(s) associated 



with resistance has been difficult, with only few exceptions [34,35]. 
Recently, microarray analysis that compares gene expression 
between susceptible- and resistant-strains unearthed candidate 
genes that confer insecticide resistance in different insect species 
[36-38]. Although molecular diagnosis of insecticide resistance 
that targets Vssc mutations is widely reported, such a system has 
not been reported for P450 genes in Aedes mosquitoes. Since 
insecticide resistance in a mosquito population could be concur- 
rently affected by more than one mechanism, an accurate 
molecular diagnosis thus needs to consider the various possible 
mechanisms. Metabolism-mediated pyrethroid resistance of A. 
aegypti requires further study especially in the populations collected 
from Southeast Asia, the largest endemic area of DF. 

Reduced cuticle penetration is the least understood mechanism 
among the three. Though it may have a primary role in resistance 
[39-41], it more often acts in combination with the other 
mechanism(s). 

This study examined the mechanisms conferring pyrethroid 
resistance in an A. aegypti strain collected from Singapore and 
investigated three major mechanisms of resistance. The results are 
expected to lead to the development of more accurate resistance 
monitoring systems and contribute to the discovery of new 
insecticide target sites to overcome the challenge of insecticide 
resistant mosquitoes. 

Methods 

Ethics statement 

Mice were used as the blood source for mosquitoes. This study 
complies with the guidelines for animal experiments performed at 
National Institute of Infectious Diseases, Japan. The protocol for 
the utility of mice was approved by the Animal Use Committee of 
National Institute of Infectious Diseases (registration numbers: 
209052, 210044, 211031 and 112029). 

A. aegypti population and strains 

A laboratory-susceptible reference strain of A. aegypti, SMK, was 
obtained from Sumitomo Chemical Co., Ltd. in 2009. This strain 
was originally from USA and has been maintained in the 
laboratory for at least 20 years without exposure to insecticides. 
The pyrethroid-resistant population, SPS 0 , was collected from 
Singapore in 2009. The SPS 0 population was selected by exposure 
to permethrin for 10 generations as described below and the SP 
strain ( = SPSi 0 ) was established. Permethrin doses used for the 
selections are listed in Table 1. Larvae were fed a ground diet of 
insect food (Oriental Yeast Co., Ltd., Tokyo, Japan), and the 
adults were maintained on 10% sucrose. Females were given blood 
meals from mice. Both larvae and adults were reared at 27 ± 1°C 
with a photoperiod of 16:8 (L:D) h. 

Chemicals 

Three technical grade insecticides were used: permethrin 
(94.1%, Sumitomo Chemical Co., Ltd., Osaka, Japan), temephos 
(94.9%, Wako Chemical Pure Industries, Ltd., Osaka, Japan), and 
pirimiphos methyl (98.0%, Wako Chemical Pure Industries, Ltd., 
Osaka, Japan). An inhibitor of P450s, piperonyl butoxide (PBO, 
98.0%, Wako Chemical Pure Industries, Ltd., Osaka, Japan), was 
used for the adult and larval bioassays. [ 14 C]-(lRS)-trans- 
permethrin (specific activity, 4.44 GBq/m mol, purity >99%) 
labeled at the phenoxy benzyl group and 4'-hydroxy-3-phenoxy- 
benzyl 3-(2,2-dichlorovinyl)-2,2-dimethylcyclopropane-carboxyl- 
ate (4'-HO-permethrin) were generously provided by Sumitomo 
Chemical Co., Ltd. (Osaka, Japan). [ 14 C]-(lRS)-trans-permethrin 
was purified by thin layer chromatography (TLC) according to the 
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Table 1. Selection of the SP strain of Aedes aegypti in the laboratory by permethrin. 



Selection No. Selected population 


Gender 


Permethrin dose (ng/mosquito) 


No. mosquitoes treated 


Mortality (%) 


1 SPS 0 G3 


0* 


10 


600 


52.0 




9 


60 


600 


62.8 


2 SPS,G1 


0* 


27 


480 


60.6 




9 


177 


480 


76.9 


3 SPS 2 G1 


0* 


44 


480 


57.1 




9 


333 


490 


85.3 


4 SPS 3 G1 


0* 


61 


513 


31.8 




9 


444 


551 


69.1 


5 SPS4GI 


cr 


111 


515 


74.0 




9 


444 


524 


50.2 


6 SPS5G2 


cr 


167 


629 


89.2 




9 


821 


482 


79.0 


7 SPS 6 G1 


cr 


167 


527 


71.3 




9 


821 


510 


78.4 


8 SPS 7 G1 


cr 


167 


583 


82.7 




9 


821 


613 


86.8 


9 SPS 8 G1 


cr 


167 


542 


65.1 




9 


821 


512 


73.6 


10 SPS 9 G1 


cr 


167 


529 


46.3 




9 


821 


547 


45.2 


The subscript numbers in the population name 


indicate the permethrin selections in the population. 


For example, SPS 5 G2 indicates the 2nd generation of SPS5 which 



was selected by permethrin for five generations. 
doi:1 0.1 371 /journal.pntd.0002948.t001 



modified method described by Shono et aL [42]. We used a 
solvent system of toluene /carbon tetrachloride (1:1) instead of 
benzene/carbon tetrachloride (1:1) in order to separate trans- 
permethrin from other degraded compounds. 

Permethrin selections and adult bioassay 

Three- to seven-day-old adult male and virgin females were 
used for the permethrin selections. In order to avoid measurement 
error caused by insect locomotive activity, we used a topical 
application method, instead of filter paper method, for pyrethroid 
selections and adult bioassays. The mosquitoes were anaesthetized 
with CO a and placed in a 200-ml Erlenmeyer glass flask. They 
were then re-anaesthetized for 3 min with 100 ju.1 diethyl ether 
absorbed onto cotton. A volume of 0.22 J_ll permethrin in acetone 
solution was dropped onto the thoracic notum of each mosquito 
with a Hamilton repeating dispenser PB-600 equipped with 10 |ll 
syringe 701SNR (Hamilton, Reno, NV, USA). The insects were 
placed in groups of 40 in paper cups sealed with a nylon mesh. 
Water-absorbed cotton was placed on the top of the mesh. 
Approximately 1000 mosquitoes were selected from each gener- 
ation (Table 1). The doses of permethrin used for the selections 
were determined by a preliminary experiment for each gender. 
Twenty-four hours after treatment, all the mosquitoes were put 
into a cage, provided with sucrose water, and maintained for 
the next generation. A bioassay to monitor permethrin 
susceptibility was conducted for each selected generation as 
described below. 

Adult bioassays were conducted by topical application using 3- 
to 7-day-old female mosquitoes and consisted of at least five doses 
of permethrin causing >0% and < 100% mortality. All bioassays 
were run at 25°C with 40 mosquitoes per dose. Mortality was 



assessed 24 h after treatment, with mosquitoes that could not stand 
on the bottom of cup counted as dead. The LD 50 values for each 
strain were calculated using log-probit mortality regression 
analysis [43]. To estimate the contribution of P450s to resistance, 
a 0.22 |il (5 M-g) aliquot of piperonyl butoxide (PBO) was placed on 
the thoracic notum of the mosquitoes 24 h prior to dosing with 
permethrin. A preliminary experiment confirmed that this PBO 
dose did not kill any of the mosquitoes. The resistance ratios were 
calculated based on the ratio of LD 5n values for the SMK and SP 
strains. 

Crossing 

To obtain Fl progeny, 100 virgin females and 100 males were 
randomly crossed and the offspring was used in the permethrin 
bioassays. BC1 progeny was also obtained by backcrossing 100 
females Fl (SPO x SMKCT) and 100 SMKff. Opposite backcross- 
ing was also conducted. Permethrin bioassays were conducted for 
both BC1 progenies as described above. 

Larval bioassays 

Larval bioassays were carried out as described elsewhere [44]. 
Twenty early fourth instar larvae were exposed to different 
concentrations of insecticides in 50 ml distilled water and 
mortality was counted after 24 h exposure at 25°C. Alcohol 
solutions of the insecticides (x200 of final cone, 250 pi) were 
added to the water. For estimating synergistic effects, the larvae 
were treated with a sublethal dose of PBO (5 |J.g/ml) in 
combination with permethrin. Mortality was assessed 24 h after 
treatment, with larvae that could not swim to the surface counted 
as dead. Three replications were run for each insecticide 
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concentration, and the LC50 values for each strain and insecticide 
were calculated using log-probit mortality regression analysis [43]. 

Genomic DNA extraction 

DNA was extracted from individual male mosquitoes using the 
REDExtract-N-Amp Tissue PGR Kit (Sigma, MO, USA). We 
used male mosquitoes to avoid gene contamination derived from 
sperm DNA in females. In A. aegypti, sex is determined by M-locus 
which is mapped on the 1 st chromosome [45] . Since Vssc gene is 
located on the 3rd chromosome [24], theoretically, genotyping of 
Vssc genes from males does not cause any gender bias. The 
extraction and tissue preparation solutions were mixed, and each 
mosquito was homogenized in a 200-ul PCR tube for 1 min at 
1500 rpm using a shaking homogenizer (MM300, Retsch Co., 
Haan, Germany) in 62.5 ju.1 of the mixture containing a zirconia 
ball (4 mm in diameter). After homogenization, the samples were 
incubated at room temperature for 30 min, followed by incubation 
at 95°C for 3 min. A 50-jal aliquot of neutralization solution was 
then added and mixed by vortexing. The extracted mixture was 
stored at -20°G. 

Genotyping of knockdown resistance gene 

Male mosquitoes were genotyped for their Vssc alleles as 
described previously [46] . We targeted six amino acid positions to 
identify the candidate(s) for knockdown resistance: the typical kdr, 
L 1 0 1 4F, and five amino acid positions in Vssc previously identified 
from pyrethroid resistant A. aegypti (i.e., S989P, 1101 1M or V, 
V1016G or I, F1534C, and D1763Y) [16,17,19,20,22]. In this 
study we numbered the amino acid position according to the 
sequence of the most abundant splice variant of the house fly Vssc 
(GenBank accession nos. AAB47604 and AAB47605). Partial 
DNA fragments of domains II, III, and IV were amplified by PCR 
using TaKaRa Ex Taq Hot Start Version (Takara Bio, Shiga, 
Japan) and three primer sets: AaSCF20 and AaSCR21 (for 
domain II), AaSCF7 and AaSCR7 (for domain III), and A1SCF6 
and A1SCR8 (for domain IV). The primer sequences are listed in 
Table SI. The cycling conditions for PCR were as follows: initial 
denaturation at 95°C for 5 min, followed by 35 cycles of 94°C for 
30 s, 55°C for 30 s, 72°C for 1 min, and a final extension step at 
72°C for 5 min. The PCR products were then treated with illustra 
ExoStar (GE Healthcare UK Ltd., Litde Chalfont, England) to 
remove unincorporated primers and dNTPs, followed by sequenc- 
ing with the following primers: AaSCF3 (forward primer for 
domain II), AaSCR22 (reverse primer for domain II), AaSCRB 
(reverse primer for domain III), and A1SCF7 (forward primer for 
domain IV) using an ABI 3130 Genetic Analyzer (Applied 
Biosystems, Foster City, CA, USA). The sequences were 
assembled and aligned using the GENETYX software (GENE- 
TYX Corporation, Tokyo, Japan). 

In vivo metabolism 

A dose of 600 dpm (ca 0.88 ng) [ 14 C]-permethrin in 0.22 ui 
acetone was administered to the thoracic notum of female 
mosquitoes. Each treated mosquito was isolated in a scintillation 
vial (6 ml) with a water-absorbed small cotton pad. The 
mosquitoes were anaesthetized with diethyl ether at 0.75, 1.5, 3, 
6, 12, 24, and 48 h after treatment and rinsed in methanol (0.5 ml 
x 2). The methanol was then added to 4 ml of scintillation 
cocktail Ultima Gold fluid (PerkinElmer Inc., Waltham, MA, 
USA) and analyzed in a liquid scintillation counter LSC-3100II 
(Hitachi Aloka Medical, Ltd., Tokyo, Japan). To study effect of 
synergist on permethrin penetration, PBO (5 u.g/9) was pl ace d on 
the thoracic notum or thoracic sternum 1 h before permethrin 
treatment. To assess the effects of solvent on permethrin 



penetration, 0.22 |0.1 of acetone was applied to the thoracic notum 
1 h prior to permethrin treatment (negative control). Each time 
point was replicated four times, with 32 mosquitoes (four 
mosquitoes x eight time points) being used for each experiment. 

After being rinsed in methanol, each mosquito was placed into a 
2-ml safe-lock tube (Eppendorf Co., Ltd., Hamburg, Germany) 
with a zirconia ball (4 mm in diameter) and homogenized for 
1 min at 1500 rpm with 0.5 ml of scintillation cocktail Ultima 
Gold fluid (PerkinElmer Inc., Waltham, MA, USA) in a shaking 
homogenizer (MM300, Retsch Co., Haan, Germany). The 
homogenized solution was transferred to a scintillation vial and 
rinsed again with 1 ml of scintillation cocktail. An additional 2.5- 
ml of cocktail was added to the vial (total volume 4 ml) and the 
radioactivity inside the mosquito was counted using a liquid 
scintillation counter. 

To measure the rate of permethrin excretion, the holding vial 
was extracted with 4 ml of scintillation cocktail and the 
radioactivity was counted using a liquid scintillation counter. 
Each time point was replicated four times (four individual 
mosquitoes) in each experiment. The rate constants for the in 
vivo dynamics of topically applied permethrin were calculated 
using the linear one compartment model described previously 
[47,48]. 

To identify the metabolites excreted into the vial, 1 0 SP females 
had ca 17 ng (11,800 dpm) of [ 14 C]-trans-permethrin in 0.44 Jul 
acetone applied to the thoracic notum and were then placed into 
20 ml glass vials. Three replicates were conducted. Forty-eight 
hours after treatment, the mosquitoes were anaesthetized with 
diethyl ether, removed from the vial, and the residual isotope in 
the vial was extracted with methanol (3 ml x 2). The methanol 
was evaporated by N 2 gas, and aliquots (5000 dpm) were spotted 
onto silica gel plates 60 F254 (HPLC, 0.2 mm, Merck KGaA, 
Damstadt, Germany), followed by development in a solvent of 
toluene/ethyl acetate (6:1). The HPTLC plate was auto-radio- 
graphed, scanned with a Bio-Imaging Analyzer BAS2500 (Fuji 
Photo Film, Tokyo, Japan), and the signal intensities of each 
metabolite were quantified. Two-dimensional development was 
conducted for water soluble metabolites using a solvent of 
chloroform/methanol/water (65:25:4) and auto-radiographed as 
described above. Unlabeled authentic 4'-HO-permefhrin was co- 
chromatographed and identified by viewing under ultraviolet light 
at 254 nm. 

Preparation of microsomes 

Microsomes were prepared from the abdomens of 3-7-old 
female mosquitoes using a modified procedure as described 
previously [49]. Mosquitoes were anaesthetized with C0 2 and 
transferred to a glass Petri dish placed on crushed ice. Two 
hundred abdomens were separated from the thoraces of mosqui- 
toes using two sets of forceps and put into a glass container on ice 
containing 5 ml of homogenization buffer [49]. The collected 
abdomens were homogenized with glass-teflon Dounce homog- 
enizer for 20 strokes using homogenizing mixer HK-1 (Asone, 
Osaka, Japan) and filtered through a layer of nylon wool. 
Homogenates were then centrifuged at 4°C at 10,000 xg for 
15 min in an Eppendorf 5804R (Eppendorf Co., Ltd., Hamburg, 
Germany) equipped with a F-34-6-38 rotor. The 10,000 xg 
supernatant was centrifuged again at 4°C at 100,000 xg for 1 h in 
a Beckman Optima MAX-XP (Beckman Coulter Inc., Brea, CA, 
USA) equipped with a MLS-50 rotor. The microsomal pellets 
were resuspended by homogenizing in 2 ml resuspension buffer 
[49]. The supernatant of the 100,000 xg centrifugation was also 
collected and used for an in vitro metabolism study. 
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In vitro metabolism 

In vitro metabolism was studied as described previously with 
some modifications [44]. The 2-ml reaction mixture contained 
incubation buffer [0.1 M sodium phosphate buffer (pH 7.5) 
containing 1 mM EDTA, 0.1 mil DTT, and 1 mM PMSF 
dissolved in ethylene glycol monomethyl ether], microsomes or 
supernatant equivalent to 10 abdomens, 0.2 ml of 10 mM (3- 
NADPH, and 100,000 dpm (0.147 ug) of [ 14 G]-trans-permethrin 
in 10 (il ethanol. An incubation mixture without [3-NADPH 
served as the control. For the inhibition study, 1 0 (0.1 of PBO 
(20 mM) or 4'-HO-permethrin (0.5, 5.0, and 10.0 Hg/ml equiv- 
alent to x7, x70, and xl40 final concentration of permethrin, 
respectively) in ethanol was added to the mixture. The mixtures 
were incubated for 5, 30, 60, 120, and 360 min at 25°C with 
shaking, and then incubation was terminated by adding 0.2 ml of 
1 N HC1, Mowed by 1 g of (NH 4 ) 2 S0 4 . Each sample was 
extracted with diethyl ether (4 ml x 3), by vortexing for 1 min and 
then centrifuged at 4000 xg for 1 min, evaporated by N 2 gas, and 
redissolved in 100 ul methanol. Aliquots of the extracted 
compounds (5,000 dpm) were spotted on HPTLC plates and 
developed in a solvent of toluene/ethyl acetate (6:1). The HPTLC 
plate was auto-radiographed, scanned with a Bio-Imaging 
Analyzer, and the signal intensities were measured as described 
above. Each time point for each strain was replicated three times 
using different enzyme sources. All the experiments involved the 
simultaneous use of the SP and SMK strains. The metabolites 
being stuck at the spotting position of the plate after development 
were collected, extracted with methanol (4 ml x 2), evaporated, 
and then spotted onto HPTLC plates. The high polar compounds 
were developed with a solvent of chloroform/methanol/water 
(65:25:4) and auto-radiographed as described above. 

Microarrays 

The microarray used in this study was designed using the 
Agilent eArray platform (Agilent Technologies, CA, USA) and 
contained 60 mer oligo-probes representing > 15,000 A. aegypti 
Liverpool transcripts identified in the genome project (https:// 
www.vectorbase.org/ content/ aedes-aegypti-liverpooltranscriptsaa 
egll3fagz). Each probe was spotted at least two times at different 
positions on each array. A probe of each P450 gene was spotted six 
times. Microarrays in a 4 x 44 k format were constructed using 
contract manufacturing carried out by Agilent Technologies. The 
entire design of the 44 k array used in this study is available from 
the NCBI Gene Expression Omnibus (GEO) site as accession # 
GPL 17604. 

The strains were reared in parallel in order to minimize 
variation resulting from breeding conditions. For each life stage 
(larvae, adult males and females), four RNA sources were prepared 
from each of the SP and SMK strains reared in different trays and 
cages ( = four biological replicates). Total RNA was extracted 
from 10 fourth instar larvae or 10 three-day-old adults using 
ISOGEN (Nippon Gene Co., Ltd., Tokyo, Japan). Genomic DNA 
was removed by digesting the total RNA samples with DNase I 
using TURBO DNase (Life Technologies Co., Carlsbad, CA, 
USA). The quality and quantity of total RNA were assessed by 
spectrophotometry using Nanodrop ND-1000 (Thermo Fisher 
Scientific Inc., Waltham, MA, USA) and a bioanalyzer MultiNA 
(Shimadzu Co., Kyoto, Japan). The purified RNA was mixed with 
the internal control RNA supplied in a RNA Spike-In Kit (One- 
color, Agilent Technologies). Fluorescein-labeled cRNA were 
synthesized via a double-stranded cDNA intermediate using a 
Low Input Quick Amp Labeling Kit (Agilent Technologies). The 
cRNA derived from the SP and SMK strains were differentially 
labeled with cyanine-3 (Cy-3) dye-conjugating CTP (PerkinElmer 



Inc., Waltham, MA, USA). cRNA was then purified with Qiagen 
RNeasy Mini Kit (Qiagen, Venlo, Netherlands), and the overall 
efficacy of cRNA synthesis and fluorescein-labeling was measured 
using a Nanodrop ND-1000 spectrophotometer. The labeled 
cRNA were pooled and hybridized to microarray probes in a 
hybridization oven at 65°C for 17 h under rotation at 10 rpm 
using a Gene Expression Hybridization Kit (Agilent Technologies, 
Santa Clara, CA, USA). After hybridization washing, the 
fluorescent dyes were stabilized against ozone oxidization with 
the Gene Expression Wash Buffer (Agilent Technologies) and the 
microarray plate was then dried in nitrogen gas. The fluorescence 
of Cy-3 on each spot was scanned using a DNA Microarray 
Scanner G565BA (Agilent Technologies, Santa Clara, CA, USA). 
Spot identification and quantification were performed using 
Feature Extraction Software v. 7.5 (Agilent Technologies, Santa 
Clara, CA, USA) in the default setting. A linear and LOWESS 
algorithm (a combination of the linear method and traditional 
LOWESS method) was used for dye normalization. Flagged spots 
(such as saturated, low intensity, and statistical outlier) were 
ignored in the final data analysis. 

The ratio of transcription levels for each gene in the microarray 
experiment for the SP and SMK strains was called the 
"microarray ratio." A representative microarray ratio in each 
hybridization experiment was expressed as the average of all spots 
for a gene on an array, where the spot number for a gene is 
calculated as "unique probe number" x "replicating spot number 
(ideally n = 4)." The P-value for every spot was calculated by 
Agilent's Universal Error Model. The raw results of the 10 
microarray hybridizations are available from GEO (series 
accession^ GSE50069). The genes over expressed in the SP 
strain were selected using a cut-off of > 3-fold relative change in 
expression and a P<0.01. 

Sequencing analysis of P450 genes 

In order to design primers for real time quantitative PCR, a 
partial or full length sequence of the P450 genes (CTP6BB2, 
CYP6Z7, CTP6Z8, CYP9M4, CTP9M5, CTP9M6, CTP9M7) and an 
internal control gene (Ribosomal protein S3 gene, RPS3) were 
amplified and sequenced from eight male mosquitoes for each of 
the SP and SMK strains. We focused on these P450 genes 
according to the following criteria: over expressed in females (>5) 
and in males (>3) in SP relative to SMK in microarray analysis. 
CYP9M7 was also analyzed as it forms a gene cluster with 
CTP9M4, CTP9M5, and CYP9M6 within 65 kbp on the same 
supercontig (1.29) and is structurally related to these genes. The 
primer sequences and their regions are shown in Table SI and 
Figure SI. Genomic DNA was extracted from individual male 
mosquitoes as described above. The PCR was conducted using 
high fidelity polymerase (PrimeSTAR GXL DNA Polymerase; 
Takara Bio Inc., Shiga, Japan), treated with illustra ExoStar (GE 
Healthcare UK Ltd., Little Ghalfont, England) to remove 
unincorporated primers and dNTPs, followed by direct sequenc- 
ing with the primers listed in Table S 1 . Because two haplotypes of 
CYP9M6, namely CYP9M6vl and CTP9M6v2, were identified from 
the SP strain, the primer sets, 9M6F3 1 /9M6R2 1 and 9M6F23/ 
9M6R35 for CYP9M6vl, and 9M6F30/9M6R22 and 9M6F18/ 
9M6R35 for CYP9M6v2 were used to amplify each variant, 
followed by direct sequencing of the PCR products (Figure S2 and 
Table SI). 

Quantitative real-time PCR 

cDNA was synthesized using total RNA isolated for microarray 
analysis, and the QT' primer (Table SI) and reverse transcriptase 
ReverTra Ace (Toyobo, Osaka, Japan) according to the manu- 



PLOS Neglected Tropical Diseases | www.plosntds.org 



5 



June 2014 | Volume 8 | Issue 6 | e2948 



Pyrethroid Resistance in Aedes aegypti 



facturer's instructions. Real time quantitative PCR was performed 
using a PikoReal Real Time PCR System (Thermo Fisher 
Scientific Inc., Waltham, MA, USA). The PCR primers were 
designed using Primer Express software (Applied Biosystems, 
Foster City, CA, USA). Each PCR reaction of 10 ul final volume 
contained 5 ul SYBR Premix Ex Taq II (Takara Bio Inc., Shiga, 
Japan), 1 ul cDNA (equivalent to 10 ng total RNA), 0.4 ul of each 
forward and reverse primer (10 uM, Table SI), and 3.2 ul 
ddH20. The PCR reactions were performed under the following 
conditions: 95°C for 2 min, followed by 40 cycles of 95°C for 10 s, 
and 60°C for 30 s. The 2 AAGt method [50] was used to quantify 
the relative expression level of P450s, with RPS3 gene acting as the 
internal control. For each gene analyzed, serial dilutions of cDNA 
showed that the efficacy of amplification for all P450 genes and the 
RPS3 gene were >0.998. Three replicates of the PCR reactions 
were performed for each sample (technical replications) and each 
experiment was repeated four times using an independent RNA 
source (biological replications). 

Identification of gene copy number 

Gene copy number was determined by quantitative PCR using 
the same primer sets described above. Genomic DNA was 
individually extracted from eight virgin females of each strain 
using Get Pure DNA Kit-Cell, Tissue (Dojindo Molecular 
Technologies, Inc., Kumamoto, Japan) according to the manu- 
facturer's instructions. The quantity and quality of the DNA was 
assessed using a spectrophotometer Nanodrop ND- 1 000 (Thermo 
Fisher Scientific Inc., MA, USA). The DNA samples were diluted 
to 100 ng/ul and 1 ul used as a template. Data were normalized 
using RPS3 gene and the exact single gene copy number reported 
by the genome project. Three replicates of the PCR reactions were 
performed for each sample (technical replications) and each 
experiment was conducted using 8 individual DNA source 
(biological replications). 

Association between permethrin excretion and CYP9M6 
and CYP6BB2 genotypes 

One hundred virgin males and females of the SP and SMK 
strains were cross-mated. The offspring were then interbred to 
obtain F2 mosquitoes. Three-day-old F2 females were dosed with 
a sub-lethal amount of [' C] -permethrin (ca 600 dpm, 0.88 ng) as 
described above. Each treated mosquito was isolated in a 
scintillation vial and the amount of radioisotope excreted was 
quantified 24 h later. Genomic DNA was isolated from six legs of 
each mosquito as described above using a REDExtract-N-Amp 
Tissue PCR Kit. Both CTP9M6 and CTP6BB2 were genotyped. 
Genotyping of CTP9M6was carried by real time quantitative PCR 
using three sets of primer: 9M6F88/9M6R89 (common), 
9M6F95/9M6R97 (specific to SP), and 9M6F96/9M6R97 
(specific to SMK). The primer sequences are listed in Table SI. 
Genotyping of CYP6BB2 was carried out by analysis of the melting 
curve of different peaks in the post-PCR assay. PCR was 
performed using a TaKaRa Ex Taq Hot Start Version (Takara 
Bio, Shiga, Japan), 6BB2F21/6BB2R22 primers, 0.5 ul EvaGreen 
(Biotium, Inc., Hayward, CA, USA), and 1 ul genomic DNA in 
10 ul of reaction mixture. The cycling conditions for PCR were as 
follows: initial denaturation of 95°C for 1 min, followed by 50 
cycles of 95°C for 10 s and 60°C for 5 s. This reaction produced 
43 bp PCR products (Figure S3). The PCR and dissociation 
analysis were performed using a PikoReal Real Time PCR System 
(Thermo Fisher Scientific Inc., Waltham, MA, USA). Total RNA 
was isolated from the remaining bodies of the mosquitoes and 
cDNA was synthesized as described above. Real time quantitative 
PCR was performed for CYP9M6 and CTP6BB2 to quantify the 



levels of transcription of these genes, as described above. Forty-two 
females (24 females for each crossing) were used for the analysis. 

Expression of P450s, reductase, and b s in Sf9 cells 

Full-length cDNAs encoding CYP9M6 and CTP6BB2 from the 
SP strain were cloned into the pPSC8 protein expression vector 
(Wako Chemical Pure Industries, Ltd., Osaka, Japan) using unique 
restriction sites positioned in the vector regions (Xbal and PstI for 
CYP9M6 and Xbal and Kpnl for CTP6BB2). Because we identified 
two variants for CTP9M6 in the SP strain {CTP9M6vl and 
CTP9m6v2), these two genes were amplified using primers 
9M6F85/9M6R86 for CTP9M6vl and 9M6F85/9M6R87 for 
CTP9M6v2. For amplification of CYP6BB2 cDNA, 6BB2F19 and 
6BB2R20 primers were used and cloned into the pPSC8 vector 
(Table SI and Figure SI). A 200-ul aliquot of SfflOOII medium 
(Life Technologies Co., CA, USA) containing 2 ug of the vectors 
(pPSC8/CYP9M6 or pPSC8/CYP6BB2), 85 ng Linear AcNPV 
(the baculovirus Autographa californica nuclear Polyhedrosis 
Virus), 5 ul Insect Genejuice Transfecfion Reagent (Merck 
KGaA, Damstadt, Germany) was transfected to expressSF+ cells 
(1 x 10 6 cells in 25 cm flask, Wako Chemical Pure Industries, 
Ltd., Osaka, Japan). The infected cells were incubated at 28°C for 
five days, centrifuged at 3000 xg for 30 min at 4°C, and the 
supernatant used for protein expression. 

Full-length cDNA encoding NADPH cytochrome P450 reduc- 
tase (PRE) and cytochrome bs cDNAs was also amplified with the 
aegREDF9/aegREDR5 and aegb5Fl/aegb5R2 primer sets, 
respectively, and then cloned into the pFastBacl vector (Life 
Technologies Co., Carlsbad, CA, USA) using multiple cloning 
sites (StuI and Xbal for PRE and EcoRI and Xbal for b 5 ). 
Recombinant constructs (pFastBacl /PRE and pFastBacl/b5) 
were used to transform MAX Efficiency DHlOBac competent 
cells (Life Technologies Co., Carlsbad, CA, USA). Recombinant 
bacmid DNA was isolated according to the manufacturer's 
instructions. The Sf9 cells were infected with recombinant 
baculovirus in Grace's Insect Medium with Cellfectin II (Life 
Technologies Co., Carlsbad, CA, USA), according to the 
manufacturer's instructions. 

For preparation of cell lysates expressing P450, 50 ml of 
confluent cells were infected in a 125-ml flask, and incubated at 
28°C with shaking at 1 30 rpm. Twenty-four hours after infection, 
hemin (in 50% ethanol and 0.1 N NaOH) was added to the media 
to a final concentration of 2 ug/ ml, followed by further incubation 
for 48 hours. The cells were then centrifuged at 3000 xg at 4°C for 
30 min and the pellets were used for the in vitro permethrin 
metabolism study. 

Metabolism of permethrin by heterologously expressed 
CYP9M6 and CYP6BB2 

The cell pellets prepared as described above were resuspended 
in 5 ml of homogenization buffer [49] and sonicated with an 
ultrasonic processor Sonics Vibra-Cell VCX130 (Sonics & 
Materials, Inc., Newtown, CT, USA) for 2 min on ice water 
(pulse on 10 s, pulse off 10 s, with 20% amplitude). The 
homogenates were centrifuged at 10,000 xg for 15 min and the 
supernatant was then centrifuged at 100,000 xg for 1 h as 
described above for the preparation of microsomes. The 
microsomal pellets were resuspended by homogenizing in 2 ml 
of resuspension buffer [49], and the protein concentration was 
determined using Bradford's reagent (Coomassie Plus Protein 
Assay; Thermo Fisher Scientific Inc., Waltham, MA, USA). The in 
vitro permethrin metabolism study was conducted as described 
above using 7 mg of microsomal proteins. The experiment was 
replicated for tree times. 
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Results 

Establishment of a permethrin-resistant strain 

A field population of A. aegypti collected from Singapore was 
selected for permethrin resistance by subjecting a wild type 
population (SPS 0 ) to the chemical for 10 generations (Table 1). 
Initially, the resistance ratio (RR) of SPS 0 was 35-fold. The RR 
rose rapidly over the generations and the established SP strain 
( = SPS 10 ) developed 1650-fold resistance after 10 generations of 
permethrin selection (Figure 1A). The RR for each generation is 
summarized in Table 2. To determine the degree of dominance of 
the resistance [51], the SP strain was crossed with a susceptible 
reference strain (SMK) and bioassays were performed on the 
progenies. The degrees of dominance for Fl (SP Q x SMKo.) and 
Fl (SMK, x SPq.) were -0.337 and -0.383, respectively 
(Figure IB). 

The toxicity of permethrin was markedly increased in the SP 
strain by pretreatment with the P450 inhibitor piperonyl butoxide 
(PBO), with the RR decreasing from 1650- to 33-fold. The 
synergistic ratio (SR) of the SP and SMK strains was 1.7 and 83, 
respectively, indicating a significant role of P450s in permethrin 
resistance of SP (Figure 1C and Table 2). The low SR of SPS 0 , at 
only 5.3, suggested that the P450-mediated resistance mechanism 



was initiated or enhanced during the process of permethrin 
selection. We observed an approximately three-fold difference in 
the LD 50 s between SPS 0 (9.8 ng/female) and SP (30 ng/female) in 
the presence of PBO. This suggested that a mechanism, other than 
P450s also contributed to the development of high level of 
permethrin resistance in SP (Table 2). The SP strain exhibited a 
4.9-fold cross-resistance to pirimiphos methyl after permethrin 
selection. 

Although the laboratory selections were conducted during the 
adult stage, larvae of the SP strain also developed a high level of 
resistance to permethrin (Figure ID and Table 3), from RR of 
160-fold for SPS 0 to 8790-fold for SP. PBO also increased the 
toxicity of permethrin in the larval stage. The SRs of SMK, SPS 0 , 
and SP were 3.5, 8.7, and 126, respectively. The SP strain showed 
a low level of resistance to organophosphates, with the RRs of 
larvae to temephos and pirimiphos methyl being 1.5- and 4.2-fold, 
respectively (Table 3). 

Frequencies of knockdown resistance type alleles 

We genotyped five positions of Vssc, which potentially cause 
decreased sensitivity to pyrethroid insecticides. Sequencing of the 
partial DNA of Vssc identified three amino acid substitutions in the 
SPS 0 population compared with the SMK strains: S989P, 




Figure 1. Log dose-probit mortality lines of Aedes aegypti. (A) The development of permethrin resistance in a field-collected population of 
adult A. aegypti under laboratory selection conditions. SPS 0 is a field-collected population prior to permethrin selection. SPS, 0 is the strain established 
by permethrin selection for 10 generations and is designated as SP. SMK is a susceptible reference strain. (B) Dose-mortality lines of adult F1 (SP 9 x 
SMKo- or SMK, x SP CT ) and BC1 (F1 g x SMKo- or SMK, x Flo-) showing the inheritance of SP resistance. (C) Synergistic effects of PBO on permethrin 
toxicity in adults. (D) Susceptibility of larvae to permethrin and the synergistic effects of PBO. 
doi:1 0.1 371 /journal.pntd.0002948.g001 
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V1016G, and F1534C. Of these substitutions, P989 and G1016 
always appeared together, indicating that these polymorphisms are 
located on the same haplotype. All mosquitoes in the SPS 0 
population possessed either the P989+G1016 or C1534 haplo- 
types, with a frequency of 44% and 56%, respectively. After the 
first selection, the frequency of P989+G1016 increased to 79%, 
and following the second selection, the CI 534 haplotypes were no 
longer detected (Table 4). In SMK, all individuals were 
homozygous for wild-type at all five amino acid positions. The 
typical kdr mutation, L1014F, was not detected in any of the 
Singapore populations and would not be expected in this species 
because kdr in Aedes would require a 2 nucleotides change. 

Rate of permethrin reduction from cuticle 

We examined the rate of permethrin reduction from cuticle by 
measuring the disappearance of insecticide from the external 
surface of the mosquitoes (Figures 2A and B). There was no 
difference in the rate of permethrin reduction from cuticle between 
the SP and SMK strains, rather slighdy higher in SP strain, with 



penetration rate constants of 0.358 h 
tively. 



and 0.252 h , respec- 



Effects of PBO on permethrin penetration 

Some previous studies suggested that synergistic effects of 
synergists are not due to inhibition of enzymes but due to 
enhancement of penetration rate of insecticides by synergists [52— 
54] . Therefore, we investigated the effects of PBO on permethrin 
penetration (Figure 2A). Pretreatment with PBO in the SP strain 
did not enhance permethrin penetration but resulted in signifi- 
candy delayed its penetration (Figure 2B). Twelve hours after 
treatment, less than 20% of the insecticide had penetrated the 
cuticle in insects pretreated with PBO compared with 90% in 
insects without PBO treatment. A reduced rate of permethrin 
penetration was also observed even when we administered PBO to 
the thoracic sternum of mosquitoes. As shown in Figure 2B, 
pretreatment of the SP strain with acetone did not affect 
permethrin penetration. 

In vivo metabolism 

The marked synergistic effects of PBO on permethrin toxicity 
implicated P450s being involved in resistance. Therefore we 
further investigated trends of permethrin after penetration through 
the cuticle using isotopic tracer tests. It was observed that in the 
SMK strain, internal radioactivity reached a peak at approxi- 
mately 12 h after permethrin treatment, and then gradually 
decreased (Figure 2C). In contrast, in the SP strain, internal 
radioactivity reached a much earlier peak and was then rapidly 
eliminated. The percentage of the initial radioactivity that 
remained after 24 h in the SMK and SP strains was 59.5% and 
14.8%, respectively. Unlike the groups without PBO pretreatment, 
the internal dose of radioactivity gradually increased in the groups 
that received PBO pretreatment and did not decrease until about 
48 h after treatment. The internal radioactivity at 48 h with or 
without PBO in SP strain was 42.7% and 4.9%, respectively 
(Figure 2C). 

Consistent with the above findings, we found that the SP strain 
excreted radioactivity more rapidly than the SMK strain 
(Figure 2D). The percentage of the radioactivity that was excreted 
24 h after treatment in the SMK and SP strains was 32.7% and 
82.8%, respectively. This suggested that the SP strain can 
eliminate permethrin more efficiently. The rate constants of 
excretion (hr _1 ) for SMK and SP were 0.022 ± 0.002 and 0.101 
± 0.008, respectively. PBO significantly reduced permethrin 
excretion, suggesting that permethrin is excreted after being 
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Figure 2. In vivo distribution of [ 14 C]-permethrin topically applied to SP and SMK strains of Aedes aegypti. (A) Schematic of the 
procedure of in vivo assay. (B) Percentage of the applied radiolabel recovered in the external rinse. [ 14 C]-Permethrin (600 dpm; 0.88 ng) was applied 
to the thoracic notum. PBO (5 u,g) was also applied to the thoracic notum or thoracic sternum 1 h prior to permethrin application to observe its effect 
on permethrin penetration. Acetone was also applied 1 h prior to permethrin application to confirm whether the solvent had any side effects. Values 
are expressed as the means (±SE) for four mosquitoes. (C) Percentage of radioactivity detected in the mosquitoes. (D) Percentage of radioactivity 
excreted in the vials. (E) Thin-layer chromatogram of the compounds excreted by females of SP strain. Developing solvents: toluene/ethyl acetate 
(6:1 ). The presence of 4'-HO-permethrin was determined by co-chromatography with an authenticated chemical. (F) Percentage of permethrin and its 
metabolites excreted by females of SP strain. Values represent the average and bars represent the SE of the mean (n = 3). (G) Re-chromatography of 
the high-polar compounds stuck around origin. Developing solvents: chloroform/methanol/water (65:25:4). 
doi:1 0.1 371 /journal.pntd.0002948.g002 
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metabolized by P450s. The high performance thin layer chroma- 
tography (HPTLC) analysis of excrete from SP strain revealed that 
more than 85% of the excreted radioisotope consisted of high 
polar compounds and approximately 10% of the permethrin was 
excreted without being metabolized (Figures 2E and F). Further 
investigation using another solvent specific for water-soluble 
metabolites showed that the compounds consisted of a number 
of metabolites (Figure 2G). 

In vitro metabolism 

We carried out in vitro metabolism studies to determine whether 
the changes observed in vivo were related to the activity of P450s. 
Microsomes of the SP strain metabolized permethrin at a much 
higher rate than those of the SMK strain (Figures 3A and B). 
Permethrin metabolism by microsomes was apparendy due to 
P450 as it was NADPH-dependent and was inhibited by PBO 
(Figure 3A). Only a very small amount of permethrin was 
metabolized by enzymes in the 100,000 xg supernatant, suggesting 
that carboxyl esterase have only minimal involvement in the 
development of resistance. The rate constant for metabolism 
(min -1 ) of SMK and SP by microsomal enzymes was 0.0007 and 
0.0042, respectively. At the start of incubation, 4'-HO-permethrin 
was the major metabolite in SP, although its percentage did not 
increase. In contrast, the percentage of high polar metabolites 
increased over time (Figure 3B, Table S2). The metabolites located 
around the origin of the HPTLC plates were collected and further 
developed with chloroform-based solvent (Figure 3C). The high- 
polar compounds consisted of a number of metabolites and were 
consistent with those in the in vivo study. 

In order to assess the effect of 4'-HO-permethrin on permethrin 
metabolism, we performed an inhibition assay by incubating 
different doses of unlabeled 4'-HO-permethrin with [ 14 C]- 
permethrin. This showed that unlabeled 4'-HO-permethrin 
inhibited effective permethrin metabolism, possibly by acting as 
a feedback inhibitor (Figure 3D). 

Microarrays 

A microarray that contained 60- mer probes for all genes 
identified by the genome project [55] was constructed. We then 
compared the levels of transcription in the SP and SMK strains of 
the adult males and females and the fourth instar larvae. In this 
study, we focused on P450s and their related genes because our 
investigations including bioassays and in vivo and in vitro 
metabolism studies showed that this metabolic enzyme is the key 
factor in the development of resistance. The genes over expressed 
in the SP strain were selected using a cut-off of > 3-fold relative 
change in expression and a P-value<0.01 (Figure 4 and Table 5). 
Nine P450 genes (CTP9M6, CTP9M5, CYP6Z8, CYP6Z7, CYP9M4, 
CYP6BB2, CYP6F3, CYP6F2, and CYP4C50) were over expressed 
in both male and female insects. Out of nine, four (CYP9M6, 
CYP9M5, CYP6Z7, and CYP9M4) were over expressed in the larval 
sample as well. Cytochrome b 5 was also over expressed in all three 
samples. CYP9M6 was the only P450 that had a >20-fold 
transcription level with a significandy low (<0.01) P- value in all 
three samples. This gene had not been reported to be associated 
with insecticide resistance. Fifteen of the 22 P450 genes listed in 
Table 5 have previously been reported to be over expressed in one 
or more of pyrethroid-resistant strains/populations collected from 
Mexico [30,56], Cuba [57], Thailand [30], Martinique, Bora Bora 
[58], and the Grand Cayman [57]. Two of these, CYP9J26 and 
CYP9J28 have been previously shown to have the ability to 
metabolize permethrin [59]. Other than P450s, no glutathione 
transferases (GSTs) or carboxyl esterases, except two GSTs, was 
over expressed in SP (under a cut off of >3-fold, P<0.01). GSTD4 



(AAEL001054) was 11.9 (P = 1.6 x 10 6 ) and 7.9 (P = 0.0025) 
times more expressed in SP of female and male, respectively, 
compared to SMK. GSTD5 (AAEL001071) in SP females was 5.5 
(P = 0.0020) times more expressed than in females of SMK. 

Sequencing analysis of P450s 

In order to design primers for real time quantitative PCR, we 
used eight mosquitoes individually to obtain the sequence of 
partial or full length DNA for seven P450 genes {CYP6BB2, 
CYP6Z7, CYP6Z8, CYP9M4, CYP9M5, CYP9M6 and CYP9M7) 
and an internal control gene, ribosomal protein S3 (RPS3). We 
focused on these P450 genes according to the following criteria: 
overexpressed in females (>5) and in males (>3) in SP relative to 
SMK in microarray analysis. CYP9M7 was also analyzed as it 
forms a gene cluster with CYP9M4, CYP9M5, and CYP9M6 within 
65 kbp on the same supercontig (1.29) and is structurally related to 
these genes (Figure 5F). The full length cDNA and deduced amino 
acid sequences of CYP6BB2 were identical between the SP and 
SMK strains. However, on the ninth nucleotide after the stop 
codon, we observed a C1533T replacement in the SP strain 
(Figure S3). All eight SP mosquitoes possessed two CYP9M6 genes, 
designated as CYP9M6vl (accession number AB840269) and 
CYP9M6v2 (accession number AB840270) (Figures S3 and S4). 
The nucleotide and amino acid identities of these alleles were 
97.6% (1563/1602) and 98.9% (527/533), respectively. In the 
SMK strain, we found three CYP9M6 variants; two mosquitoes 
were homozygous for CYP9M6v3 (accession number AB846835) 
or CYP9M6v4 (accession number AB846836), two were homozy- 
gous for CYP9M6v5 (accession number AB846837), and four were 
heterozygous CYP9M6 variants. The CYP9M6v5 protein ap- 
peared to be incomplete as an active enzyme because there was a 
stop codon at amino acid position 473 (Figure S4). 

Quantitative real-time PCR 

We performed real-time quantitative PCR (qPCR) for seven 
P450 genes to verify the results of the microarray analyses. We 
focused these genes according to the criteria described above. 
Primers for quantitative real time PCR were designed from 
sequences that were consensus among eight SP and SMK 
mosquitoes. Five of these genes, CYP6BB2, CYP6Z7, CYP9M4, 
CYP9M5, and CYP9M6 were over expressed during the larval 
stage and in adult males and females (Figures 5A— C) in SP relative 
to SMK. As shown in Figure 5D there was a good correlation 
between the results of the microarrays and qPCR (R 2 = 0.824). 
This confirmed that the results of the microarray were accurate. 
CYP9M6 showed the largest relative change in the seven genes and 
was over expressed 22.9-fold in males, 38.9-fold in females, and 
28.2-fold in larvae of the SP strain relative to SMK. The relative 
expression level of CYP6BB2 to RPS3 was also higher (close to 
CYP9M6], especially in adult females (Figure 5A), with the gene 
reported to be over expressed in five other pyrethroid-resistant 
strains (Table 5). 

Identification of the gene copy number 

The gene copy numbers of the seven P450s in the SP and SMK 
strains were compared by qPCR (Figure 5E). RPS3 gene was used 
to normalize the data. While the majority of the P450 genes in the 
SMK strain were likely to be a single copy, some genes in SP were 
amplified significantly. The average copy number of CYP6Z7, 
CYP9M4, CYP9M5, and CYP9M6 in SP were 3.1-, 4.3-, 4.8-, and 
4.6-fold more than SMK, respectively. This clearly shows that the 
over expression of these P450 genes, at least in part, was due to 
gene amplification. 
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Figure 3. In vitro [ 14 C]-permethrin metabolism studies. (A) Thin-layer chromatogram of in vitro permethrin metabolism in adult Aedes aegypti 
with or without NADPH or PBO. Microsomes (100,000 x g pellet) and the soluble fraction (100,000 x g supernatant) were used as the enzyme source. 
The presence of 4'-HO-permethrin was determined by co-chromatography with an authenticated chemical. (B) Metabolites of permethrin in adult A. 
aegypti at various incubation times. For detailed data, see Table S2. (C) Thin-layer chromatogram of high-polar compounds (SP) located at the origin 
of the first HPTLC. Developing solvents: chloroform/methanol/water (65:25:4). (D) Inhibition of permethrin metabolism by 4'-HO-permethrin. 
Microsomes of the SP strain were incubated for 30 min with [ 14 C]-permethrin, NADPH, and x7- x70-, or x140-fold higher doses of 4'-HO-permethrin 
than permethrin. 

doi:1 0.1 371 /journal.pntd.0002948.g003 



Association between permethrin excretion and CYP9M6 
and CYP6BB2 

In order to determine whether CYP9M6 and CYP6BB2 have a 
role in permethrin metabolism, F2 progeny (interbred population 
of Fl (SMK x SP)) were treated with a sub-lethal dose of [ 14 C]- 
permethrin, and the excretion rate was measured 24 h after 



treatment in individuals. The genotypes and transcription levels of 
CYP9M6 and CTP6BB2 were determined for each mosquito to 
examine the association between these parameters and the rate of 
excretion (Figures 6A and D). The transcription levels for CTP9M6 
between RR (homozygote of SP allele) and SS (homozygote of 
SMK allele) were significantly different (P<0.0001), with a 45-fold 
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Figure 4. Microarray screening of cytochrome P450 and b 5 genes differentially expressed in SP and SMK strains. Differential 
transcription of cytochrome P450 and b 5 genes was individually investigated in adult female (A, B), adult male (C, D), and fourth instar larvae (E, F) of 
Aedes aegypti. Figures 4A-C show volcano plots of the relative changes (log 2 , x-axis) and statistical significance (-log 10 (P-value), y-axis). Figures 4D-F 
show scatter plots of the relative expression level of each P450 and b 5 gene. Genes with relative expression >3-fold and a P-value<0.01 are 
highlighted. 

doi:1 0.1 371 /journal.pntd.0002948.g004 



difference between the two genotypes (Figure 6B). The rate of 
permethrin excretion was 69% for RR and 45% for SS 
(Figure 6C). On the other hand, although the transcription level 
of CYP6BB2 was significandy different between RR and SS (P = 
0.0298), the relative change (RR/SS) was only 1.6-fold (Figure 6E). 
The rate of permethrin excretion in RR (7 1 %) was significant than 
that in SS (49%) (P<0.0001, Figure 6F). These results suggested 
that both CYP9M6 and CYP6BB2 were associated with rapid 
excretion of permethrin, although it was likely that the metabolism 
rate of permethrin was different between the two isozymes. 

In addition to CYP6BB2 and CYP9M6, the transcription levels 
of other four P450 genes (CYP6Z7, GYP6Z8, CYP9M4, and 
CYP9M5) were quantified in 48 F2 progeny in order to determine 
their association with the rate of permethrin excretion (Figure 7). 
We observed a relatively high correlation between the rate of 
permethrin excretion and mRNA levels for CYP6BB2 and 
CYP9M6, with correlation coefficients (R 2 ) of 0.454 and 0.485, 
respectively (Figures 7A and F). The R for CYP9M5 was also 
high (0.537, Figure 7E). Of the six P450 genes, CYP6Z7 had the 
lowest R 2 (0.236, Figure 7B). Forty-eight individuals were ranked 
from 1 to 48 according to the level of expression of each P450 
gene. The ranks of each individual against multiple P450 genes 
were combined to standardize the level of mRNA and we 
evaluated the association between permethrin excretion rate and 
expression level of multiple P450 genes. The standardized 
transcription levels of CYP9M6+CYP6BB2 had a considerably 
higher R 2 (0.569) than individual values alone (Figure 7G). This 
strongly suggests that both these P450s contribute to permethrin 
excretion. The combination of CYP9M6, CYP6BB2, and 
CYP9M5 further increased the R 2 to 0.589 (Figure 7H), whereas 
the combination of all six P450s decreased the R 2 to 0.499 
(Figure 71). These findings imply that not all of these P450 were 
involved in permethrin excretion. 

Metabolism of permethrin by heterologously expressed 
CYP9M6 and CYP6BB2 

In order to determine whether CYP9M6 and CYP6BB2 are 
capable of permethrin metabolism, these P450s were co-expressed 
with A. aegypti cytochrome P450 reductase and b 5 in Sf9 cells using 
a baculovirus. For CTP9M6, two genes, CTP9M6vl and 
CTP9M6v2 were expressed, as sequencing analysis showed that 
all individuals in the SP strain possessed both types of genes. 
Microsomes were prepared from cells infected with baculovirus, 
followed by an in vitro metabolism study using [ 14 C] -permethrin as 
the substrate. Both CYP9M6vl and CYP9M6v2 demonstrated 
relatively low but consistent permethrin metabolism, whereas 
CYP6BB2 exhibited strong metabolic activity for permethrin 
(Figure 8, Table S3). The compound, 4'-HO-permethrin, detected 
in the microsome metabolism studies, was also confirmed as a 
major metabolite for the three P450 samples (CYP9M6vl, 
CYP9M6v2, and CYP6BB2). This suggested that these three 
P450s can detoxify permethrin. In addition, a large amount of 
polar metabolites (E) were detected by in vivo and in vitro studies on 
all three samples. Unique metabolites (metabolites A and B), which 
were almost undetectable in the in vitro microsomal study were also 
detected in samples of CYP9M6vl and CYP9M6v2 (Figure 8, 
Table S3). 



Discussion 

We selected a Singapore colony of A. aegypti (SPSo) in the 
laboratory and established a SP ( = SPS 10 ) strain, which developed 
an extremely high level of resistance to permethrin. Using 
bioassays with a synergist and in vivo and in vitro studies, we 
confirmed that P450s play very important roles in the develop- 
ment of resistance. The larvae of the SP strain also developed a 
high level of resistance despite the selections being conducted on 
adult mosquitoes. This indicates that common mechanism(s) 
causes resistance in both developmental stages. This finding is in 
contrast to those reported in the JPal-per strain (Culex quinque- 
fasciatus), with these insects being selected by permethrin in the 
larval stage and not showing a high level of resistance in the adult 
stage [60] . The JPal-per strain has a knockdown resistance gene 
allele (L1014F) as well plus high expression levels of CTP9M10, 
known to confer resistance [61-63], being approximately 250-fold 
higher than in the larvae of susceptible strain, and showing limited 
expression in the adult stage [37]. In the SP strain, CTP9M6, 
CTP6BB2, and CYP9J26 were over expressed in both the larval 
and adult stages and probably these were involved in the larval 
resistance as well. 

The amino acid substitutions of Vssc, V1016G and F1534C, 
have been shown to strongly correlate with pyrethroid resistance in 
A. aegypti [17,23-27]. The Singapore population (SPS 0 ) possessed 
both G1016 and C1534, with a gene frequency of 44% and 56%, 
respectively. However, all C1534 haplotypes were disappeared 
after only two selections by permethrin. This clearly indicates that 
Vssc with a G1016 mutation has a lower sensitivity to permethrin 
than channel with a CI 534 mutation. This finding is consistent 
with that reported by a recent neurophysiological study [27]. The 
results of our bioassays further supported this phenomenon, as we 
showed 3-5-fold differences in resistance ratios between SP and 
SPS 0 insects when they were treated permethrin with PBO in the 
adult (Table 2) and also larval stages (Table 3). Changes of Vssc 
genotype is also associated with development of resistance in the 
larval stage in SP. The reason for relatively high frequency of 
CI 534 mutation in this country is still unclear. Further studies, 
including the examination of fitness costs between G1016 and 
CI 534, will help us better understand factors affect the equilib- 
rium state of these two haplotypes in nature. 

In the bioassays, PBO exhibited a marked synergistic effect on 
the toxicity of permethrin in the SP strain (Table 2 and 3). On the 
other hand, several research groups stated that synergists may 
enhance the toxicity of insecticides because they accelerate the 
penetration of these agents [52-54]. We therefore examined the 
effects of PBO on permethrin penetration and showed that it did 
not enhance permethrin penetration, but instead suppressed the 
rate of its penetration (Figure 2). This finding is consistent with 
other reports on the German cockroach [64] and house fly [65- 
67]. We also found that permethrin penetration was markedly 
suppressed by PBO treatment, even administered on the thoracic 
sternum of the insect followed by permethrin treatment on the 
thoracic notum. This result implies that PBO was not acting as a 
physical barrier to permethrin penetration but rather another 
mechanism. We speculate that high concentration of PBO inside 
the mosquito body suppressed permethrin passing through the 
cuticle by the effect of a concentration gradient of chemicals. 
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1.29 on the 2nd chromosome. The error bars represent standard errors of three (A, B, and C) and eight (E) biological replicates. 
doi:1 0.1 371 /journal.pntd.0002948.g005 



Therefore, caution must be exercised when we use a synergist in a 
bioassay to avoid underestimating the contribution of metabolic 
enzymes to the resistance. 

4'-HO-permethrin has been well documented as a primary 
metabolite of permethrin in various insects including the house 
fly, southern house mosquito, cockroach, and cabbage looper 



[42,44,68]. We also showed that this compound was a major 
metabolite of permethrin in the SP strain in vitro and comprised 
approximately 40% of the metabolites after an incubation of 
5 min (Figure 3B). However, the rate of 4'-HO-permethrin to the 
total metabolites decreased over time and comprised approxi- 
mately 10% after 6 h. On the other hand, the rate of high polar 



PLOS Neglected Tropical Diseases | www.plosntds.org 



17 



June 2014 | Volume 8 | Issue 6 | e2948 



Pyrethroid Resistance in Aedes aegypti 




TJ 

CD 

— -i 

3 

CD 



CD 
X 

o 

— 1 
CD 

5' 

Z3 



0.0001 0.001 0.01 0.1 

Relative level of CYP9M6 mRNA 



D 




o 



o ° 



-466- 



80 



60 



40 



20 



0.001 0.01 0.1 1 

Relative level of CYP6BB2 mRNA 



TJ 

CD 



CD 
X 

o 



o 




RR SR SS 

Genotype of CYP9M6 




RR SR SS 

Genotype of CYP6BB2 



^ 100 




100 



RR SR SS 

Genotype of CYP9M6 



c 

o 



o 

X 

0 



CO 

E 

q3 

Q_ 
H— 

o 

£ 
ffl 
or 




RR SR SS 

Genotype of CYP6BB2 



Figure 6. Genotype-gene expression-phenotype associations. The association between CYP9M6/CYP6BB2 genotypes and mRNA expression 
and permethrin excretion rate was investigated. The female F2 populations of SP and SMK strains were dosed with [ 14 C]-permethrin, and the amount 
of permethrin excreted was individually quantified 24 h after treatment. The genotype of each mosquito was identified by the genomic DNA isolated 
from six legs. mRNA level was also quantified with RNA isolated from the body part of each mosquito. For more information, see Materials and 
Methods. (A, D) Plot showing the association between the level of mRNA and permethrin excretion in CYP9M6 (A) and CYP6BB2 (D). (B, E) Box plot 
showing the association between mRNA level and CYP9M6 (B) and CYP6BB2 (E). (C, F) Box plot showing the association between permethrin excretion 
rate and genotype of CYP9M6 (C) and CYP6BB2 (F). P-values for comparison of the RR and SS genotypes were calculated using Dunnett's test. 
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Figure 7. Associations between the gene expression level of six P450s and permethrin excretion rate. The female F2 populations of SP 
and SMK strains were dosed with [ 14 C]-permethrin and the amount of permethrin excreted by each insect was quantified 24 h after treatment. The 
level of mRNA was quantified using RNA isolated from each mosquito. For more information, see Materials and Methods. (A) CYP6BB2, (B) CYP6Z7, (C) 
CYP6Z8, (D) CYP9M4, (E) CYP9M5, and (F) CYP9M6. Forty-eight individuals were ranked from 1 to 48 according to the level of expression of each P450 
gene. The ranks against multiple P450 genes were combined to standardize the level of mRNA and evaluate the association between permethrin 
excretion rate and expression level of multiple P450 genes. (G) Association between permethrin excretion and standardized level of mRNA for 
CYP6BB2 and CYP9M6. (H) Association between permethrin excretion rate and standardized level of mRNA for CYP6BB2, CYP9M6 and CYP9M5. (I) 
Association between permethrin excretion rate and standardized level of mRNA for CYP6BB2, CYP6Z7, CYP6Z8, CYP9M4, CYP9M5, and CYP9M6. 
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compounds increased over time and made up 86% of the total 
metabolites after 6 h. Further, 4'-HO-permethrin was barely 
detectable in the metabolites excreted by mosquitoes in vivo 
(Figures 2E and F). This suggested that 4'-HO-permethrin was 
further metabolized to secondary metabolites that were more 
water soluble and excre table (Figure 9). We also found that 
4'-HO-permethrin inhibited permethrin metabolism in vitro 
(Figure 3D), indicating the presence of a negative feedback 



regulation (Figure 9). Insects need to convert hydroxyl-permethrin 
to other secondary metabolites in order to avert this negative 
feedback, resulting in insufficient excretion. Secondary develop- 
ment in the HPTLC analysis showed that the high polar 
metabolite consisted of a number of various compounds 
(Figures 2G and 3C). These compounds were suspected of being 
conjugates of 4'-HO-permethrin with glucosides or/and with 
various amino acids including glycine, glutamic acid, glutamine, 
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Figure 8. Thin layer chromatogram of [ 14 C]-permethrin metab- 
olites by heterologously expressed CYP9M6v1, CYP9M6v2, 
and CYP6BB2. A-E denote unidentified metabolites. Developing 
solvents: toluene/ethyl acetate (6:1). The presence of 4'-HO-permethrin 
was determined by co-chromatography with an authenticated chem- 
ical. 

doi:1 0.1 371 /journal.pntd.0002948.g008 

and serine, as reported by Shono et al. [68] . It remains unknown 
what enzymes play the role in the secondary metabolism of 
pyrethroid in insects, and it is possible that an obscure mechanism 
of insecticide resistance may exist. 

Association of gene amplification with insecticide resistance is 
relatively well documented [69]. In organophosphorus-resistant 
aphids and Culex mosquitoes, genes of carboxyl esterases are 
amplified 80-250 times and the capability to detoxify insecticides 
is increased [70-72]. Recent studies demonstrated that gene 
amplification is also associated with over expression of P450 genes 
in insects [40,62,63,73,74]. In A. aegypti, over expression of 
CYP9J29 in the CUBA-DELTA and CAYMAN strains was 
shown to be due to gene amplification [57]. In the current study, 
we found that CYP9M6 had the capability to metabolize 
permethrin and was over expressed in the SP strain partially due 
to gene amplification. Real time quantitative PCR revealed that 
this gene was amplified approximately 4-fold in the SP strain 
compared with the SMK strain. We identified two CTP9M6 genes, 
CYP9M6vl and CTP9M6v2, which were probably raised by the 
gene amplification event in the SP strain. In particular, two other 
similar genes, CTP9M4 and CTP9M5, located within the same 
cluster on chromosome two [75] were also amplified to a similar 
extent (Figure 5E). This suggested that this amplification may have 
occurred coincidendy. 

Our microarray analysis revealed that 22 P450 genes were over 
expressed > 3-fold in all samples (males, females, or larvae) of the 
resistant SP strain compared with the susceptible SMK strain 
(Table 5). It is noteworthy that 15 of these genes have been 
reported previously to be overexpressed in one or more resistant 
populations or strains collected from different regions of the world 
[28,30,56,57]. In this study, we found that CYP9M6 and 
CYP6BB2 were responsible for resistance. In addition, CYP9J26 
and CYP9J28, which have the ability to metabolize permethrin 
[59,76], were also over expressed in the SP strain. It remains to be 




4'HO permethrin 



1 Conjugation with various 
compounds (sugar? amino acids?) 

High polar, soluble, more 
excretable compounds 

I 
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Figure 9. Proposed model showing the metabolic pathway of 
permethrin. 

doi:10.1371/joumal.pntd.0002948.g009 

determined as to how many P450s contribute to the resistance in 
SP mosquitoes. However, it is not likely that all genes over 
expressed in the SP strain are involved in the development of 
resistance, as it is generally known that transcription of genes 
forming clusters are occasionally coordinately controlled by the 
same regulatory factor [77]. In association analysis using F2 
progeny, the expression level of CYP9M5 and CYP9M6 showed a 
high correlation coefficient (R~ = 0.90), whereas no correlation 
was observed between the expression levels of CYP9M6 and 
CTP6BB2 (R 2 = 0.10, data not shown). This suggested that 
CTP9M5 and CYP9M6 are controlled partially by a common 
mechanism. A further example of this mechanism is the finding 
that administration of phenobarbital to insects induces a multiple 
number of P450 genes [78,79]. Therefore, over expressed genes do 
not always provide an advantage for insect survival. Heterologous 
expression and confirmation of the metabolic activity of each P450 
isoform is necessary to gain a greater understanding of these 
processes. 

To the best of our knowledge, this is the first report showing a 
correlation between the expression level of each P450 and the rate 
of insecticide metabolism of individual insects. This analysis may 
provide a good tool for evaluating the rate of contribution of each 
P450 isoform to the development of resistance. In this analysis, 
both CYP9M6 and CYP6BB2 showed moderate correlation 
between the levels of transcription and the individual rate of 
permethrin metabolism. Furthermore, although no correlation 
was observed between the expression level of CYP9M6 and 
CYP6BB2 (R = 0.10), there was a prominent relationship 
between the rate of permethrin metabolism and standardized- 
combined expression level of CYP9M6 and CYP6BB2 (Figure 7G). 
These results strongly suggest that these two P450s are auxiliary 
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parts of effective permethrin excretion. Although CTP9M6 showed 
the greatest relative change in gene expression in the comparison 
between the SP and susceptible strains, we still hesitate to conclude 
that this enzyme plays a major role in permethrin metabolism. 
First, homozygotes of CYP6BB2 had a rate of permethrin excretion 
as high as that observed in homozygotes of CTP9M6 in Figures 6C 
and F. Second, in the heterologous expression study, under our 
experimental conditions, CYP6BB2 showed a much higher 
activity of permethrin metabolism than CYP9M6. It is therefore 
possible that CYP9M6 may compensate for its low metabolic 
activity by its high mRNA level. Although both CYP9M4 and 
CYP9M5 were also over expressed in the SP strain, we did not 
investigate the permethrin metabolic activity of these P450s 
because their relative expression level was below 1/10* of that of 
CYP9M6. However, it is possible they may confer resistance 
depending on their unit activity of permethrin metabolism. 
Further work will be needed to clarify if these P450 isoforms are 
involved in the resistance. 

One of the goals of the study was to establish a more accurate 
molecular diagnosis method for monitoring field populations of A. 
aegypti by identifying metabolic enzymes that confer resistance. 
Findings of CYP9M6 and CYP6BB2 are good progress for this 
purpose, however, we also found that the high level of resistance in 
the SP strain was not due to one major mechanism but was 
actually the consequence of multiple P450 isozymes (at least 4) and 
reduced sensitivity of Vssc. This implies that we may need to 
recognize the reality that the development of molecular diagnoses 
targeting metabolic enzymes will be more complicated and 
challenging than that of targeting only Vssc [80] . We are currendy 
focusing on elucidating the contribution degree of these P450s on 
pyrethroid resistance in the field populations of A. aegypti collected 
from different regions. 
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